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Abstract 

We formulate SU(3) Hamiltonian lattice QCD in terms of the pla- 
quette variables and determine the relevant subspaces of the Hilbert 
space for the vacuum wave functional and its approximations. We 
analyze the one- and two-plaquette problems. 

In an earlier paper || we formulated the fully gauge-fixed SU(iV) lattice 
gauge Hamiltonian and studied the one-plaquette problem for several values 
of N. This paper is restricted to the case of SU(3) lattice gauge theory. For 
this case I discuss ways to go beyond the one-plaquette approximation. 

The Hamiltonian approach to lattice QCD, although more complicated, 
has advantages over the Euclidean formulation. || [| In general, one can 
say that the Hamiltonian approach allows one to split the problem in parts 
and tackle each part separately, thereby generating more understanding and 
restricting the eventual computational effort. This should be formulated 
consistently as restrictions on the full physical Hilbert space to a proper 
subspace, or likewise, either by formal arguments, or from calculations. One 
important step in this process is the gauge fixing, such that all the degrees 
of freedom are physical, and an excitation spectrum, determined from an 
ab-initio many-body approach, such as RPA, is actually physical, and does 
not contain oscillations in the gauge parameter. The first and foremost part 
of the full QCD problem is the determination of the vacuum wave functional. 
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The dimensionless QCD lattice Hamiltonian can be split into three parts: 



" ' «x,0> 



where the sum over a runs over all plaquettes, and the sum over < a, (3 > 
runs over all nearest-neighbour plaquette pairs. The dimensionless energy e 
is related to the energy E for a given lattice spacing a and oupling constant 
9- 
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and the coupling constant A can be expressed in terms of the gauge-field 
coupling constant g: 

If we ignore the nearest-neighbour interaction, the Hamiltonian is local and 
the solution is the product wave functional of one-plaquette wave functions. 
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It is important to realize that although a special unitary 3x3 matrix U is 
parametrized by eight angles, only two are relevant for the one-plaquette 
problem. Since the potential, or magnetic, term h M || only depends on the 
eigenvalues x, y and z: 



Tr[U + rf] 





X 






' x~ l 


Tr 




y 

z 


+ Tr 


y- 1 

z- 1 


x + y ^ 


1 

■z + ~ 

X 


1 

+ - + 

y 


1 

z 



(4) 



of the plaquette variable U with the relation xyz = 1 it is sufficient to know 
x = expicj) x and y = expiry. In the case of several plaquettes, we can 
simultaneously diagonalize multi-plaquette states with up to five different 
plaquettes, only for correlations of six and more plaquettes simultaneous 
diagonalization might not be possible in special cases. 

For the ground state at zero temperature there are no excitations in the 
transverse directions because the Hamiltonian does not couple to these states, 
so the wave functional is constant for these angles.Q So for our analysis of 

lr The normalization is not a problem since the group is compact. Correspondingly, the 
maximal- Abelian gauge is not fully Abelian because of correlations of sixth and higher 
order spoil the diagonalization in three spatial dimensions. 
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the vacuum, we can restrict ourselves to the variables x and y, however, for 
symmetry reasons, we use also z, keeping in mind that z = (xy)~ x . 

So in terms of the eigenvalues x, y, and z of the plaquette variable U the 
one-plaquette Hamiltonian is: 



h = 2A + 4S - ^(T 2 + 6T)-^(x + y + z + - + - + - 
3 2 \ x y z j 
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where we define 
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The operators A and T 2 are diagonal for each monomial x p y q z l , and give 



the respective eigenvalues, p 



I and (p + q + I) . The terms S and 



T are due to the curvature of the group manifold given by the Jacobian 
(x — y)(y — z)(z — x). Only S mixes the different monomials, however, 
perserves the order p + q + /. It is easy to see that the electric operator 
is consistent with the SU(3) group, since for any eigenstate <f>(x,y,z) one 
can show that (xyz) r <f)(x, y, z) is an eigenstate with the same eigenvalue, 
which then can be mapped on (j>(x,y,z). In the strong-coupling limit A = 
0, the solutions are the eigenstates of the electric operator, which form an 
orthonormal basis on the group manifold. Generally, we do not have to 
perform any integration over the group, as the eigenstates of the electric 
operator guarantees orthogonality of the basis states. Actually, all we need 
for an eigenvalue equation is a unique labeling of states, orthogonal or not, 
an invariant subspaces, on which the action of the Hamiltonian is bijective, 
such that a truncation to a finite subspace is possible. A similarity transform 
can make the Hamiltonian hermitian if necessary. For the one-plaquette 
problem these states are orthogonal and characterized by a label (n, m) where 
n > m > 0: 
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with eigenvalues |( 



pair (n, m) 



n" +m 2 —nm — 3) . The eigenvalues are degenerate for each 
and (n, n — m), which corresponds to the hermitian conjugate 
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(n, m) 


X(n, m ) = E Cijkmijk, i + j + k = n + m-3, i>j>k 
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Table 1: The eigenstates (77,777) expressed in monomials, for example: 
m 3,2,2 — x 3 y 2 z 2 + y 3 x 2 z 2 + z 3 x 2 y 2 and 7772,2,2 — x 2 y 2 z 2 . These expressions 
can be reduced, using the relation xyz = 1, however, for most purposes it 
useful to keep the original form. The last column contains the corresponding 
eigenvalues. 



solutions. In principle, for the vacuum wave functional, the basis can be 
restricted to symmetric combinations only: ^(|n, m) + |n, 77 — m)). There 
are several equivalent representations of these states. For example, they 
can be expressed as polynomials of trace variables: t\ = x + y + z and 
t 2 = x 2 + y 2 + z 2 , or in terms of loop variables: = (x z +y l + z l ){x^ +yi + . 
In the mappings between these representations one makes extensive use of 
the relation xyz = 1. In order to derive results efficiently we need to go to 
and fro between representations. We choose the representation in terms of 
symmetric polynomials of the eigenvalues because of its clarity and versatility. 
I will return to this representation to check signs and factors in the formulae. 

In Table [l] we express the eigenstates up to 77 is six in terms of monomials. 
The, generally integer, coefficients of the mapping between representations 
can be deduced from group theoretical methods using Young tableaux, || 
but in this paper I will not refer to any group-theoretical result. A general 
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Figure 1: The occupation numbers of the lowest basis states for the one- 
plaquette ground state. Due to the symmetry between covariant and con- 
travariant representations, the occupation numbers for the states (n, m) and 
(n, n — m) are identical. 



formula for an electric eigenstate is given by: 

(x,y,z\n,m) = J2 J2 ^-^ n+m ~ 2_i ~V ' 2 1 ~ l ~ k y ' 2 * " +k 

= X] C i,j,k m i,j,k , (8) 

% > j > k 
i+j+k=n+m— 3 

which follows from the evaluation of the determinants of Eq. [F[ We used the 
result, which will be useful later for the evaluation of the 5 operator, 

J'' 1 -' 1 ' = (xyr^ ' l ]T ^MI-^V • (9) 

x y k=o 

In all this, o is the sign function, with o"o = 0. 

If we use this basis up to a certain number^] n we can solve the one- 
plaquette problem for arbitrary A = g~ A . Away from the strong coupling 
limit a large basis is required for convergence, but with 100-200 states a 



The Hamiltonian restricted to this subspace is hermitian. 
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good convergence is obtained for A < 50. In Figure [I] we show the occupation 
number for the lowest states as a function of the coupling constant. Clearly, a 
typical expansion around the strong coupling limit with, say, 5-10 states will 
break down quickly, as many electric eigenstates contribute to the ground 
state. In Figure § we show the lower end of the spectrum. 

Therefore, in order to go beyond the one-plaquette approximation, it is 
not only necessary to construct multiple plaquette states, but also restrict 
the physical Hilbert space to make any calculation feasible. A straight two- 
plaquette problem requires already thousands of states for a reasonable con- 
vergence, and there will be no easy way to extend that to larger lattices. We 
will return to this problem later. 

From the one-plaquette problem, we are guided to write the Hamiltonian: 

H = y £h a + £ E*E% , (10) 

a <a,/3> 

where h a is the Hamiltonian for the plaquette a, which we solved earlier. Al- 
though the ground-state and the low-lying states of the one-plaquette prob- 
lem contain many contributing basis states, it would not be too hard to con- 
struct an effective Hamiltonian with only the eigenstates of the lower end of 
the one-plaquette spectrum. The real question is how to construct the spatial 
excitations; the important multiple-plaquette correlations, which then in its 
turn are used to construct the short-range interaction in the effective Hamil- 
tonian. In a way, part of the physics ends up in the coupling-constant and 
lattice-spacing dependence of the parameters of this effective Hamiltonian. 

Since the operator E^E^ = VL a p + V a /3, contains a diagonal part fl a/ 3 and 
a part V a p between plaquette a and plaquette /3, which brings us out of the 
one-plaquette product wave functional. We can use it as generator for spatial 
excitations. For a state \(n a ,m a ))\(np,mp)) the multiple application of the 
operator V a/ 3 can generate up to mm{n a , rip} — 2 spatial excitations. In order 
to see this come about we look at a basis state in terms of its trace variables. 
One term may be written as: 

W 3l32 — 'ji = £ Ui x i-Ui 2 i^- ■ 'Uiji^ (11) 

where {jij'2.73 • • ■ ji} is an permutation of {123 ••■/}. In the case of the iden- 
tity: {J1J2J3 • • • ji} = {123 • • • 1} it is a product of single traces: Wj = Tr[U] 1 
and for {jij2j3 • • • ji} = {^12 •••(/ — 1)} it yields Wj = Tr[U 1 ]. So in the end it 
is not the permutation {,71,72,73 • • • ji} that gives a unique representation, but 
the cocycles in the permutation; the number of traces and each exponent. 
This can be labeled by exponents as the partition of I. For the special group 
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Figure 2: The low-lying energy spectrum for the one-plaquette ground state. 

SU(iV) the representation in terms of trace variables is not unique. Because 
xyz = 1 every term can be expressed in terms of two traces Tr[[/*]Tr[£P] since 
a product of three traces contains terms like x l y^z k from which a common 
factor [xyz) mm ^ 1, ^ k ^ can be removed. This has been the greatest headache for 
people studying SU(3) Hamiltonian lattice gauge theory^. It mixes also 
terms of different order. As soon as one includes more complicated states, in 
terms of trace variables, the abundance of redundant states is much larger 
that the number of unique states. There does not seem to be an easy way to 
determine the unique and complete set in terms of trace variables. 

However, we will work from the other way around; we have the electric 
eigenstates of the SU(3) group on one plaquette, without redundancies. We 
construct the corresponding trace variables from that, and we can now study 
the effect of the correlation operator V a p. We start best from the SU(3) 
identity: 

1 8 1 1 

^ Yl KjKl = 2 6il6k 3 ~~ gM« ■ (12) 

The application of the electric operator corresponds to the insertion of these 
A- matrices at the appropriate places. So acting on the right, for example, on 
the trace above Eq. [11] at places, r and s, we find that 

1 1 

{jl •••Ji}- > ^{ir" jr-ljsjr+1 ■ ■ ■ js-ljrjs+1 ' ' ' jl} ~ ' ' ' jl} ■ (13) 

2 o 
3 For SU(2) this problem is much easier. 



As simple as that. 

We have the electric operator on each plaquette, which is a combination 
of the permutation operator and the identity operator for each basis state, 
and we have the inter-plaquette operator V a /3 that permutes lines between the 
two plaquettes. So in the two-plaquette problem we distinguish six operators: 
P a and P/3 the permutation operators of lines of plaquette a and plaquette 
P respectively, Q a and the diagonal part of the electric operator on each 
plaquette, V a p the operator that permutes lines of plaquette a and plaquette 
(3 , and Q, a p the diagonal part of the inter-plaquette operator. 

Now we can make a number of important observations. Firstly, every 
eigenstate on a plaquette of the electric operator can be expressed in terms 
a finite number of normal- ordered permutation operators P without common 
indices on a given number of lines I: 

n+m—3 [5] 

(U\n,m)= :«•:■■■> , (14) 

1=0 k=0 

where the normal-ordered permutation operator : P fc : can be expressed in 
terms of permutations P(i,j) of elements i and j, and the inter-plaquette 
operator V accordingly: 

:P": = E P{iiM)P{hM)---P{i2k-iM) , (15) 

ii 7^2 i z ---i z iik 

■V k ap : = E P(iiJi)P(i2,j2)---P(ikJk) , (16) 

where the indices ii, - • • ,ik belong to the plaquette variables on the plaquette 
a, and ji,---,j fc to the plaquette f3. Terms like P(12)P(23)wi 23 = Tr[f/ 3 ] 
can be reduced to w + |P(12)wi 2 3 — ^123 = w o + |P^i23 — |^i23, where 
w = Tr[l]. 

Secondly, the inter-plaquette operator can be normal ordered, without 
common indices, and truncates at the minimal number of lines on one of 
the plaquettes. The normal ordering generates additional inter-plaquette 
permutations: 

v aP v a ,wt ] ,A% = ■■ vie ■■ ™^A% (1 7 ) 

+ pqw\..'. p w\..'. q . 

Lastly, the magnetic term on each plaquette can be recovered from rela- 
tions like: 

Tr[U] : P k : = L-J^ : p* : (18) 
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W]:P* = <Ll^Z»+a : p*« : ro ,.., ( , +1)(H2 , (19) 

So every two-plaquette basis state can be represented by five labels: 

({U a ,U fi }\k,p,l,q,r) =: 7^ : : P£ : : PJ : ufiWJ , (20) 

where p, q, and r are restricted by k and Z, such that p < fc/2, g < Z/2, 
and r < min{fc,/}. Linear combinations of : P p : Wi...k span the same 
space as linear combinations of (n, m) with n + m — 3 = k. An important 
corollary of this result is that different one-plaquette electric eigenstates do 
not mix due to the inter-plaquette action, since only one eigenstate with a 
given eigenvalue lies within the space spanned by : P p : Wi...^. 

In terms of electric eigenstates on each plaquette there would be a re- 
dundance, since, for example, : V^p '■ ({U a Up}\(3i 2) a (4, 1)^) = due to 
cancellations among the trace variables, however, due to mixing of states 
with different n there is no real redundancy, if the commutation relations 
are worked out. Similarly, a straightforward truncation of the I < L and 
k < K would violated hermiticity, and it would lead to complex eigenvalues, 
since the symmetry between the covariant (n, m) and the contravariant rep- 
resentation (n, n — m) is broken by this truncation. It is easy to mistake the 
non-permuting part of Eq. for the diagonal part of inter-plaquette part of 
the electric operator, however, some care is required here since the basis in 
trace variables is not orthogonal. 

We will illustrate most of the concepts above in a simple example. Let's 
consider one-plaquette diagrams with two lines. Starting from Tr[?7] 2 = w\2 
there are two diagrams: w\2 and Pwi2- The electric operator in this case is: 

E 2 = , (21) 



3/ 3 

where the last term is the diagonal part from Tr[U] J2 a Tr[A a A a t/], ^ and the 
first part is the specific case of Eq. |TJ where P = 2P(1,2) = P(l,2) + P(2, 1). 
The eigenvalue equation is given by: 

* ( ^ ) . (22) 




2 



PWl2 



which leads to the eigenstates (3, 2) and (4, 1) and the eigenvalues: A( 3j2 ) 
8/3 and A(4 i i) = 20/3 as denoted in Table |]. 



4 This equals (N 2 — 1)/N for N = 3 times the number of matrices U in the product 
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If we now look at the two-plaquette problem, we see that the one-plaquette 
electic operator commutes with the inter-plaquette part: 



[P a ,Pp] = [P ,V aP ] = [P a ,V a/3 ] = 



(23) 



therefore we can work independently in the four separate subspaces of {0(3 <Pu 

{0(3 2)) 0(f i)}- The two-plaquette electric operator has a simple action on 
each of these states 



E 1 



(24) 



where the factor | is because the permuting part of the electric operator on 
the common link is: \(P a + P/3 + TP a fi). On all the other links there is a 
single iP/3 or \P a . For the mixed states we get the eigenvalue equation: 




^(3,2)^(4,1) 
2 ' a/3V ) (3,2)V ) (4,l) 



x | ^(3,2)^(4,1) 

2' a/3</ ) (3,2) ( / ) (4,l) 



(25) 



while for the symmetric states the electric operator yields: 




( <f> ( f ] \ 

l<r> A a P) 

\ 4 • ' a/3 ■ VA J 



( 



A 



(a/3) 
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(a/3) 



da/3) 
>A 



(26) 



where A = 38/3 and the positive sign for (4,1) states of each plaquette, 
i.e., 4>^ = 0(4i)0(fi), an d A = 14/3 and the negative sign for (3,2). The 
two largest of the three eigenvalues of the two-plaquette electric eigenstates 
for (3,2)0,(3,2)0 are 14/3 and 17/3, which correspond with the eigenval- 
ues of coupled (3,l) a (3,l)0 states. A representation in terms of symmet- 
ric polynomials in x a ,y a and z a , and xp,yp and zp reveals indeed that 
these are hermitian conjugate solutions. The lowest eigenvalue 8/3 of the 
(3,2)0(3,2)0 triplet is degenerate with the lowest single-plaquette excitation 
like (3,2)o(2, 1)0, however, this corresponds to a rebundant solution, as the 
three states of Eq. ^B] are linearly dependent for A = (3, 2), contrary to the 
case of A = (4, 1). 

This procedure can be carried out to construct any eigenstate of the 
electric operator without generating many redundant states. Furthermore, 
it does not rely on the diagonalization of plaquette variable U in the one- 
plaquette problem, so it can be applied hierarchically up to any order in a 
multi-plaquette Hilbert space. However, in the future we hope to map this 
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procedure completely onto a basis of symmetric polynomials which should 
remove the last redundancies. 

Eventually, as we increase the spatial extent of the states, we can restrict 
ourselves to the lower energy states of the subsystems from which these 
states are build. Thereby we restrict the computational effort, which, as 
mentioned before, is the main problem of an ab-initio Hamiltonian lattice 
QCD calculation. Finally, I gratefully acknowledge discussions with Niels 
Walet during the early stages of this work. 
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